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TITLE 



TYPE 



Post Multiplication of a Matrix A or Its Transpose 

rrr 

A "by a Vector k. 

Closed with one program parameter «. Used with N - 8. 

The link'Iss 

A0j?F 



NUMBER OF 
TEMPORARY 
PARAMETER 



26 xF 
where A is the location beginning at which non-zero 

elements of the matrix A are stored and x is the first 

T 
word of this routine. Set A = K for Ak; A = 5 for A k, 

WORDS kj 

STORAGE 0,1 

S3, S4, S5- During the input of this routine locations 

3, k and 5 must hold the following quantities s 

n is the numb ere of rows in the 

matrix A. 

Location (b + i), (i = 0, <».«,, q-l), 

contain the q components of the 

vector k. 

Locations (c + l)(i = 0, o«o, q-l) 



3) 00 F OQr-nF 
k) 00 F 00 bF 



5) 00 F 00 cF 



contain the q components of the 

,T, 



ACCURACY 
DURATION 



vector product Ak or A k, 
A maximum of 9 decimal places 



r 

z 

i=l 



(1.53 c. + o98) ms for Ak, 



£ m 

L (1-59 c. + 080) ms for Ak, 

1=1 1 



- 2 



where 



c. = number of non-zero elements in row r, 
r = number of rows of A. 

DESCRIPTION This routine will postmultiply any matrix A or its 

T 
transpose A by a cblumn vector k. It is designed to handle a matrix where 

only the non-zero elements are stored. These non-zero elements are assumed 

to be stored consecutively.. Each non°zero element utilizes the last ten 

bits to indicate the row and column in which it occurred o The first non-zero 

-30 
element of each row has one in the 2 ^ -th location! and the other non-zero 

-30 
elements of the same row have zero in the 2 th -location. The remaining 

-31 -39 
nine bits 2 s <>»<>, 2 indicate the column index » The last non-zero 

-30 
element of A must be followed by a word containing 2 o Routine H - 8 

will input and store the non-zero elements of a matrix properly modified <» 

LIMITATI01S This routine will not handle any matrix with more 

than 511 columns . The elements a . . of the matrix 

A and the elements k. of the vector k must be scaled 

J 

so that 



E lA. . k. I < 1 in the case of Ak. 



and 



V I T | T 

L A.„ . k. < 1 in the case of A k« 

i=l ° 

This will be true if 

ii/ 2 

|k. J < x/n with x < n 



and 



|a„ . I < x/m in the case of Ak. ot 

i T i / T 

a„ . < x/m' in the case of A k v 
1 ij ' ' 9 



-3- 



where 

n a number of components of the vector k 

m = max (r^ r^, . .., r , ..-) 

where r . is the number of non-zero elements in 

the i th row. 

m« = max (c^, c 2 , . .., c , ...) 

where c is the number of non-zero elements in 

the i th column. 

It is advisable to scale down by powers of 10 so that 

one can count back merely by shifting the decimal 

point . 

&0TES l) If one wishes to perform only Ak type multi- 

plications, he may delete the last 20 words of 
this routine, and thus reduce it in length to 
25 words . 

2) If one wis.es to use this routine repeatedly 
he should clear out the locations (c + i) between 
entries. Otherwise, the components of the previous 
multiplications iLk, will have been added to the 
respective components of the present multiplication 
" : A 9 k with the final result being A.k + A k instead 
of A 2 k 2 * 
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42 12L 






1 


36 25L 




T 
to A k 




10 20F 




_ Set addresses for finding 


2 


42 17L 




J elements of A 




42 20L 






3 


4l LF 








Fl IF 




Set counter to -1 


4 


40 IF 




. 




L4 5F 




Set address for 


5 


42 15L 
00 20F 




storing Ak components 


6 


46 15L 
22 17L 






7 


40 F 
50 F 
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F5 15L 








40 15L 




New Row 
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L5 11L 




reset address es,» 




l4 15L 




- increase counter 


10 


46 15L 




and test for end 




26 11L 




waste 


n 


F5 IF 
40 IF 






12 


LO 3F 
32 ( )F 






13 


S5 F 
10 24L 






14 


26 20L 
7 J F 






15 


L4 f 
40 F 






16 


F5 17L 
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40 17L 
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LO 
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22 


i4l 


00 
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00 
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IF 


LO 
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Form Ak 



Test for new row 



(1023 +512) S4 
Constants 



Set addresses for finding 

T 
elements of A 



Set counter 
= -1 



Set for address of 
1st component of k 



New Rcw s reset 
addresses,, increase 
counter^ and test 
for end 
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LO 24L 
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38 


42 42L 
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Form A k 
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S5 F 
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